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Abstract 

■ The effect of unitary noise on the performance of Grover's quantum search algorithm is studied. 

o ■ 

. This type of noise may resuh from tiny fluctuations and drift in the parameters of the (quantum) 

components performing the computation. The resulting operations are still unitary, but not pre- 
cisely those assumed in the design of the algorithm. Here we focus on the effect of such noise in 

o ■ 

CN ■ the Hadamard gate W, which is an essential component in each iteration of the quantum search 



> 
(N 



process. To this end W is replaced by a noisy Hadamard gate U. The parameters of U at each 
iteration are taken from an arbitrary probability distribution (e.g. Gaussian distribution) and are 



' characterized by their statistical moments around the parameters of W. For simplicity we assume 

O '■ 

. that the noise is unbiased and isotropic, namely all noise variables in the parametrization we use 

o , 

, have zero average and the same standard deviation e. The noise terms at different calls to U are 

assumed to be uncorrelated. For a search space of size N = 2^ (where n is the number of qubits 



used to span this space) it is found that as long as e < 0(n "^N 4), the algorithm maintains 



Q^, significant efficiency, while above this noise level its operation is hampered completely. It is also 



found that below this noise threshold, when the search fails, it is likely to provide a state that 
■ differs from the marked state by only a few bits. This feature can be used to search for the marked 

state by a classical post-processing, even if the quantum search has failed, thus improving the 
success rate of the search process. 

PACS numbers: PACS: 03.67.Lx, 89.70.+C 
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I. INTRODUCTION 



The discovery of quantum algorithms that can solve computational problems faster than 
any known classical algorithm stimulated much interest in quantum information science. 
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ie_known algorithms include Shor's factoring algorithm Grover's search algorithm 

as well as algorithms for the simulation of physical systems. One of the most serious 
obstacles that should be dealt with in the way to construct a quantum computer on which 
such algorithms can be implemented is the problem of decoherence |5|. This is the effect 
of the interaction between the quantum system that stores and manipulates the quantum 
information and the environment, that spoils the coherence of the quantum states. The 
effect of decoherence on Grover's search algorithm was recently studied using perturbation 
theory p. Recent progress in quantum error correcting codes lid . as well 

as in decoherence free sub-spaces 13, llE M, provide an effective way to keep 



the quantum states coherent and enable the implementation of useful quantum algorithms. 
However, a drawback of these approaches is that they involve some redundancy in the 
encoding of the logical quantum state, thus requiring to maintain a larger number of quantum 
bits in a coherent state. 

The performance of a quantum computer may also be affected by unitary noise 
Such noise may result from tiny fluctuations and drift in the properties of the implemented 
quantum gates. These fluctuations or drifts may add stochastic perturbation elements to 
the Hamiltonian that describes the quantum gates that generate the unitary operations in 
the algorithm. The perturbated Hamiltonian is Hermitian as well, therefore the resulting 
operations are still unitary, although not precisely the ones assumed in the design of the 
algorithm. Since the implementation of a useful quantum algorithm requires a large number 
of one and two qubit gates, it is possible that even a tiny noise in each operation would 
accumulate to a considerable effect that may hamper the operation of the quantum computer. 

In this paper we analyze the effect of unitary noise on Grover's quantum search algorithm. 
To this end we replace the Hadamard gate W hj a noisy (but still unitary) Hadamard gate U. 
The parameters of U at each iteration are taken from an arbitrary probability distribution 
(e.g. Gaussian distribution) and are characterized by their statistical moments around those 
ofW. 

In order to simplify the calculations we assume an unbiased noise (i.e. noise with mean 
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) or alternately we refer to a known bias that is shifted away in every iteration of the 
Hadamard gate. We assume that the noise is isotropic, namely the standard deviation e is 
the same for all the noise variables in the parameters to be defined later. This assumption is 
made for simplicity, and removing it does not change the main qualitative features observed 
in this paper. 

We show that the noise reduces the success probability Pq of the algorithm (the probability 
to measure the marked state, the one we are looking for, at the end of the algorithm). 
Consider a search problem with N elements where N — 2'^ {n denotes the number of qubits 
in the register) . We find that as long as the standard deviation e of the noise satisfies 



the algorithm maintains significant efficiency, while above this noise level, its operation is 
hampered completely. 

We have analyzed the flow of probability out of the marked state which is caused by the 
noise. In the effective region of the algorithm (i.e. e < eo ), we find that a considerably large 
amount of probability diffuses from the marked state to its near neighbors, whose indices 
differ only in a few bits from the marked index. This result is derived analytically and 
verified numerically. 

The "diffusive" flow from the marked state to its near neighbors can be used to enhance 
the efficiency of the algorithm. To this end we execute the Grover quantum search, and 
measure the state of the register. If it is not the marked state wc classically test all its 
neighboring states, namely those that differ from it by only a few bits. The use of hybrid 
(quantum and classical) search strategies in case of unbiased and isotropic unitary noise 
in the Hadamard operations accomplishes a great reduction in the average searching time 
in comparison to the ordinary quantum search procedure that is re-execution of Grover's 
algorithm over and over again until the marked state is found. Moreover, hybrid strategies 
enable an efficient quantum search under noise levels for which the quantum search alone 
fails. 

The paper is organized as follows: In section II we introduce the Grover quantum search 
algorithm. In section III we analyze the effect of unitary noisy Hadamard gates on the 
Grover quantum search. In section IV we present and discuss numerical results that verify 



e<eo = 1.4 
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the analytical predictions and extend the scope of discussion to noise's limits for which the 
analytical approximations are not valid anymore. In section V we examine the use of hybrid 
search strategies which improves the performance of the Grover quantum search with noisy 
Hadamard gates. Some additional details of calculations are presented in the appendix. 



II. GROVER'S SEARCH ALGORITHM 



Let D be a search space containing N elements. We assume, for convenience, that = 2", 
where n is an integer. In this way, we may represent the elements of D using an ra-qubit 
register containing their indices, i = 0, . . . , N — 1. We assume that a single marked element 
is the solution to the search problem. The distinction between the marked and unmarked 
elements can be expressed by a suitable function, f : D ^ {0,1}, such that / = 1 for the 
marked element, and / = for all other elements. 

Suppose we wish to search the space D to find the marked element, namely the element 
for which / = 1. To solve this problem on a classical computer one needs to evaluate / for 
each element, one by one, until the marked state is found. Thus, 0{N) evaluations of / 
are required on a classical computer. However, if we allow the function / to be evaluated 
coherently, there exists a sequence of unitary operations which can locate the marked element 
using only 0(-\/iV) queries of / . This sequence of unitary operations is called Grover's 
quantum search algorithm. 

To describe the operation of the quantum search algorithm we first introduce a register, 
\x) = \xn-i, Xn-2, ■ ■ ■ , Xi,Xo) of Ti qubits, and an ancilla qubit, |g), to be used in the compu- 
tation. The bar script denotes a computational string of n bits, i.e. x = Xn~2, ■ ■ ■ ,xi,xo 
where the fc'th bit x^ is either one or zero, so that \x) is a computational basis state vector. 
We also introduce a quantum oracle, a unitary operator Of which functions as a black box 
with the ability to recognize the solution to the search problem. We shall use the hat script 
for operators, namely an operator O is denoted as O. (For more details on how an oracle 
may be constructed, see Chapter 6 of Ref. Q].) The oracle performs the following unitary 
operation on computational basis states of the register and of the ancilla \q): 

df\x)\q) = \x)\q(Bf{x)), (2) 
where © denotes addition modulo 2. This definition may be uniquely extended, via hnearity, 
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to all states of the register and ancilla. 

The oracle recognizes the marked state in the sense that if x is the marked element of 
the search space, f{x) = 1, the oracle flips the ancilla qubit from |0) to |1) and vice versa, 
while for unmarked states the ancilla is unchanged. In Grover's algorithm the ancilla qubit 
is initially set to the state 

k) = ^(|o)-|i)). (3) 

It is easy to verify that, with this choice, the action of the oracle is: 



V2 J \ V2 

Thus, the only effect of the oracle is to apply a phase rotation of tt radians if x is the marked 

state, and no phase change if x is unmarked. Since the state of the ancilla does not change, 

it is conventional to omit it, and write the action of the oracle as 

= (5) 
Grover's search algorithm may be summarized as follows: 

1. Initialize the qubit register to |0) = |0,0, . . . ,0) and the ancilla to |0). Then apply 

. (oi\ 

the gate wx on the ancilla qubit where x — I is the NOT gate and w — 

[ I is the Hadamard gate (The matrices are written with respect to the 

computational basis |0), |1)). The resulting state is: 

,o,o,...,o)(MH.)^ 

2. Grover Iterations: Repeat the following operation T times: 

(a) Apply the Hadamard gate on each qubit in the register. 

(b) Apply the oracle, which has the effect of rotating the marked state by a phase of 
TT radians. Since the ancilla is always in the state (|0) — \l))/^/2 the effect of this 
operation may be described by a unitary operator acting only on the register. 
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= Ylx(~'^)^^^^\-'^) ~ ^ ~ 2\m){m\. (j denotes the 2" x 2" identity operator 
and \m) denotes the marked state, which is the one we are searching for.) 

(c) Apply the Hadamard gate on each qubit in the register. 

(d) Rotate the |0,0, ...,0) state of the register by a phase of tt. This rotation is 
similar to 2(b), except for the fact that here it is performed on a known state. It 
takes the form Jq = -|0)(0| + ^^^^ \x){x\ = / - 2|0)(0|. 

3. Apply the Hadamard gate on each qubit in the register and then measure the register 
in the computational basis. 

We still need to specify the number of iterations, T. As subsequent Grover iterations are 
applied, the amplitude of the marked state gradually increases, while the amplitudes of the 
unmarked states decrease. There exists an optimal number, T, of iterations at which the 
amplitude of the marked state reaches a maximum value, and thus the probability that the 
measurement yields the marked state is maximal. Let us denote this probability by Pq- It 
has been shown [J, |2^ that the optimal time T is: 



T < 



(7) 



/N 

4 

where [x] denotes the integer value of x. Moreover, it was shown that Grover's algorithm 
is optimal in the sense that it is as efficient as theoretically possible . For the initial state 
given in step 1 above, t he p robability to obtain a marked state at the optimal time T is 



= i-oii/VN) \i^m |. For other initial states Pq may be bounded away from unity 
22 1 . Further generalizations to the original Grover's algorithm have also been suggested in 

"n 

Mm- 

Following Grover's discoveries, a variety of applications were developed, in which the 



algorithm is used in the solution of other problems 
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H H 0, Q, Q 3, 31, y, Q; 

Experimental implementations were also constructed using a nuclear magnetic resonance 
(NMR) quantum computer [^^Q] as well as on an optical device 0]. 



III. ANALYSIS OF THE EFFECTS OF UNITARY NOISE 

Assuming a single marked state, the original Grover quantum search algorithm using 
an n-qubit register with N = 2^ computational basis states, could be represented by the 
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operation: 

\g{T))^Go{T)\0) = WQ^\0). (8) 

\g{T)) is the quantum state after T iterations which is accomphshed by executing the Grover 
operator Gq{T) on the initial state |0) = |0 . . . 00). When the time T is optimal, the operator 
Go{T) in Eq. (jHl) amplifies the absolute value of the amplitude of the marked state to 1. 
The Grover operator Go{T) consists of T executions of Grover's iteration Q followed by a 
single Hadamard transform W. The Grover iteration is defined as the composite operation: 

Q = -foWLw (9) 

where: 

n-l n-1 I ^ ^ 

w = (^Wk = 6^^\ I (10) 



k=0 k=0 ^ 



1 -1 



k 



is the n-qubit Hadamard transform {® denotes a tensor product and (S)fcIo ~ ^n-i ® 
. . . ® wq), Jq = / — 2|0) (0| is a selective rotation by vr of the state |0) and Im = I — 2\m) {m\ 
is a selective rotation by tt of the marked state \m) (J denotes the identity operator). The 

matrix notation used is in the computational basis representation, i.e., |0)^ = I I and 







k 



|1) 



k 





1 



k 

Grover's iterations consist of the Hadamard transforms Wk, k = 0, . . . n — 1 which are 
one qubit gates, and the selective inversion operators that involve two qubit gates. Here 
we consider the effect of unitary noise in the Hadamard gates Wk, k = 0, . . .n — 1 on the 
performance of Grover's algorithm. To this end we define the noisy Grover iteration at time 
t as: 

Qt = -ioVjmUt. (11) 

Here Ut and Vt are the direct products of one qubit unitary operators, which are built of noisy 
executions of one qubit Hadamard gates. We assume that there are no temporal correlations 
in the Hadamard executions and that operations on different qubits are uncorrelated as well. 



7 



Therefore, the operators can be written as: 



ra-l 



k=0 
n-1 

T4 = (S?) e^^'^W^fe = e^^W (12) 



k=0 

where dk{t) and Ck{t) are one qubit stochastic operators, acting on the /c'th qubit in the 
register at time t. They are Hermitian and generate the unitary noise in the one qubit 
Hadamard gates. Any one qubit Hermitian operator can be expanded in the basis of the 
Pauh operators &ik, &2k and a^k and the identity operator ijofc, acting on the /c'th qubit. 
Therefore: 

3 

^k{t) = ^<^nk{t)crtj,k 
n=o 

3 

Ck{t) = ^lnkit)^nk (13) 

where: 

/^io\ ^ ( ^ ^\ ^ ( ^ ( ^ ^ 

V'l. vv. \'°). v°-^* 

(14) 

are the identity and Pauh matrices in the computational basis representation and Q;^fe(i) 
and ^ij,k{t) (a* = 0, 1, 2, 3 and k — 1, . . . , n) are real stochastic variables. Since all the ct's 
and 7's are produced by the same physical hardware with no correlations, their statistical 
properties can be expressed by their moments: 

{l,k{t)) ^ (7^) (15) 

and 

{af,k{t)auk'{t')) = 5tt'5kk'{5Q^f_i5ay) + {a^){au) 

{lu.k{t)lvk'{t')) = ^tt'Skk'^ln^lu) + {'yi^)M (16) 

where: 

= 7m - (7m) (17) 



and 5u' as well as 5kk' are Kronecker's delta functions. The noise is thus characterized by 
the eight real stochastic variables: and 7^ (yU = 0, 1, 2, 3). 

Using a series expansion of the exponentials in Eq. (jl2j) and the identity of operators: 

{A®B){C®D) =AC®Bb 

we obtain: 



n-l 



n-l 



(18) 



k=0 



The operators At and Ct are built as sums of n one qubit operators. Therefore, the only 
non-zero elements in their 2" x 2" matrix representation in the computational basis are 
those having indices that differ in no more than one bit. Denoting two computational basis 



states: |x) = Xn-2, ■ ■ ■ xq) and 



IVn-i, yn-2, • • • 2/0) {xk = 0, 1 and Uk = 0,1 where 



k = 0, 



1), the following matrix elements can be calculated using Eqs. (fT ^ - (fT^ and 



(x\A, 



EVoi^okit) + {-l)y^a,kit)} if ||x - y\ 

aifc(t) + ii-iy''a2kit) if I |x - y| I = 1 (x^ = -^yk), 
if I Ix — |/| I > 1 







(19) 



and 



EVoilokit) + i-iy^lskit)} if ||x = , 

{x\Ct\y)= I 7u.(t)+^(-l)^^72fc(t) if ||x-y|| = l ix, = ^yk), (20) 

if I |x — y| I > 1. 

The norm | |x — ?/| | is the Hamming distance that counts the number of bits that are different 
in X and y. Similarly, one can write: 



n-l 



k=0 

n-l 

wCtW = Y,^Mt)wk 



(21) 



k=0 



SO that: 



{x\WAtW\y) 



if I |x — y\ 







«3fc(i) - i(-l)^'=«2fc(t) if ||x - yll = 1 {xk = -^yk) 
if I |x — yl I > 1 



(22) 



9 



and 



' n-l 



{x\WCtW\y) = < 



EVoiiokit) + i-iy'iikit)} if = , 

73,(t)-2(-l)J"=72fe(t) ii\\x-y\\ = lixk = ^yk), (23) 
if I |x — 1/1 1 > 1 



are elements of 2" x 2" sparse matrices as well. 

The noisy Hadamard transforms lit and Vt defined in Eq. (I12j) can now be written 
explicitly into the noisy Grover's iteration for any time t: 

Qt = -he'^'WLWe'^' 

= he'^'UQe'^' = e'^°^'^'Qe'^\ (24) 

and the entire search process after T iterations is then given by: 

G{T) = Ut+iQt-.-Qi 

= (H/e^^^+i)(e^^o^^^oge^^^)(e*-^"°^^-i-^~« . . . e'^^){e'^''^'^''Qe'^'). (25) 

A Taylor series expansion now yields 

e'^^^^e'^^'^^^'^ ^ (/ + ^At+l - \a\^^ + ...)(/ + ^kCtk - \hClh + • • •) 

^ / + lEf^ - \e? + . . . (26) 

where 

Ef^ = it+i + ioCjo 

Ef~^ = i^^i + 2At+JoCth + hClh (27) 

are the first and the second order deviations, caused by the noise, so that: 

G{T) ^ W{i + ^E?) - ^4') + . . .)Q{i + lE^^l, - \eP_, + ...)... Q(/ + - + . . .) 
= Go(T) + Gi(T) + G'2(T) + .... (28) 

where we define Co = 0. The operator 

Go(T) = WQ^ (29) 
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is the original Grover quantum search operator (without noise), while Gi{T) are the com- 
ponents of perturbation of order I (l — 1,2...). Specifically, the leading terms of the 
perturbation satisfy: 

T 

G^{T)^^W{Y,{Q^-'M'''Q')} (30) 
t=o 

and 

G,{T) = -lw{f2iQ^-'EPQ') + 2 E '"^{Q^-'Ei'^Q'-'MI^Q')}. (31) 
t=o t=i t'=0 

Those additional perturbation components consist of summations over time indices t of 
operator multiplications. These multiplications include the noise generators At and Ct where 
the number of their appearances determines the order of the perturbation. They also include 
non-stochastic operators: powers of the original Grover iteration Q and selective inversions 
around the initial state of zeros Iq. In order to understand the behavior of the operator G{T), 
let us first focus on the non-stochastic operators which appear in those multiplications. 
For any arbitrary state 1-0) of the n-qubit register 

= -{i -2\^{^\)W{i -2\fh){fh\)W\'il)) 
= -|V^) + 2{(0|t/;) - 2{Q\W\m){m\W\il])}\Q) + 2{m\W\ 

In particular: 

g|0) = {1 - 4|(m|W^|0)f }|0) + 2{m\W\Q)W\m) 

and 

QW\fh) = -2{Q\W\fh)pj + W\fh). 

The operator Q acts as a linear transformation within a 3-dimensional vector space that is 
spanned by the three vectors: |0) and W^|m). In case that the state vector is linearly 
independent on the vectors |0) and H^|m), the following representation of vectors basis can 



)W\m) 



(32) 

(33) 
(34) 



be performed: 





VV 



|5> 



W\m) 





VV 



(35) 



Note that the these vectors are not represented in the computational basis. However, 
they are linearly independent and hence can be treated as a standard basis in a simple 
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3-dimensional vector space in which an inner product is not defined. Namely, the above rep- 
resentation does no^ imply orthogonality of the vectors l'^), |0) and iy|m), and the following 
matrix calculations are exact and involve no approximations at all. This representation does 
not change eigenvalues and eigenvectors from those obtained for an orthonormal basis. 

Now, denoting (0|H^|?ti) = (m| W^|0) = the operator Q, which is a regular linear trans- 
formation acting on a simple 3-dimensional vector space, is represented by the 3-dimensional 
matrix: 

^-10 



12 



4^ 1 



(36) 



where: 



qi = Him 



[m\W\'ip)} 



q2 = 2{fh\W\'ilj). (37) 

Powers of the operator Q do not exceed from the 3-dimensional sub-space spanned by \iJj), 
|0) and H^|m). Diagonalizing the matrix Q we obtain 

^Ao ^ 

Qd= A+ (38) 
^0 A_ ^ 

where the eigenvalues are Aq = — 1 and X± — e^"^, and cu is given by: 



cos a; 



1 - 



(39) 



The diagonalizing matrix whose columns consist of the three eigenvectors takes the form: 

f 1 ^ 
1 1 



X = 



Vl 

V2 e+^^ e-''^ 



(40) 



where 



and 



2 



1-e 



(41) 



V2 



1 - 



N 



- {m\W\i^) 



1-^ 



(42) 
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We now obtain an explicit expression for powers of Q 



/ 



(-1)' 

vi[{-i)' -k{t)]+v2m k{t) -i{t) 

yv2[{-iy - m{t)] - vil{t) l{t) m{t) j 



(43) 



where: 



m 
m 

m{t) 



sin(ci;t — ip) 

sin((y9) 
sin(co't) 



sin(v9) 
sin(ct;t + ip) 
sin((y9) 



(44) 



The expressions above can be simphfied in the hmit of large N . In this limit the frequency 
uj given by Eq. (jHTIjl can be approximated by 



N 



N2 



(45) 



Under this approximation Eq. (|4H) becomes 



namely 



As a result: 



and 



TT 



k{t) m{t) = cos(u;t) + O 



(46) 
(47) 



l{t) = sin(cjt) + O 

Moreover, as we shall see later, the superposition lip) appears in our first order calculation 
in only two forms. Either lip) = \y) where \y) is a computational basis state perpendicular 
to |0) or l^/j) = W\y) where \y) is a computational basis state perpendicular to the marked 
state |m). In the first case we obtain 

^{m\W\y) ^ ( 1 



V2 



-{m\W\y) 
1 - i 



o 
o 



N 



(48) 
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In the second case we find that 



V2 



-(0|^l^) 

^(01^1^) 
1-^ 



o 



(49) 



Thus in both cases the matrix of Eq. ()43p takes the form: 

/ 



{-If 











0(-^) COS (cjt) — sin (cjt) 
\'-'^Vn^ sin(^t) cos(^t) J 



(50) 



This means that acts simultaneously, up to 0{^^), as a 2-dimensional rotator in the 
sub-space spanned by the state vectors |0) and W\rh), as well as selective phase invertor 
(according to the parity of the power t) of any state vector {tp) that is independent on |0) 
and W^|m). 

The selective inversion around |0) (denoted as Jq) also preserves the 3-dimensional simple 
vector space spanned by {ip), For any superposition lip) 



io\^P) = {i-2\0){Om 



-2(0|V^)|0), 



(51) 



where {0\ip) = for the case of {ip) = \y) 7^ |0) and {0\ip) = for the case of \ip) = W\y) 
{\y) 7^ \fn)), where \y) denotes an arbitrary computational basis state. Particularly we can 
write: 

/o|0) = -|0) (52) 

and 



ioW\m) = W\m) - -^\0) 



N 

Hence in both cases the selective inversion can be represented by: 



(53) 



v 



1 

1 



(54) 



The effect of the stochastic operators At and Ct, given by Eqs. (fTUI) - on a compu- 
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tational basis state \y) is given by: 

n-l 

At\y) = {y\At\y)\y) + Y,{yik\At\y)\yik) 

k=0 

n-l 

AtW\y) = + 

A:=0 

n-l 

Ct\y) = {y\Ct\y)\y) + Y,{yik\cmyik) 

k=0 

CtW\y) = + (55) 

k=0 

Here {||/ifc)}o<fc<n-i the n basis vectors, that are different from the basis vector \y) in 
a single qubit only, denoted as the fc'th qubit in the register. For example, in a 5 qubit 
system, where \y) = |17) = |10001), {\yik)}o<k<4 computational basis states 

which includes |10000), |10011), |10101), |11001) and |00000). Any operation of either 
At or Ct on basis state vectors of the form \y) or increases the dimension of the 

relevant vectors space by n additional independent directions determined by {||/ifc)}o<fc<n-i 
and {W\yik)}o< 

k<n-i respectively. 

Grover's search algorithm is performed by executing the operator G{T) on the initial 
state of zeros |0). Ignoring the effect of the unitary noise in the Hadamard transforms (i.e. 
considering G{T) ^ Go{T) = WQ'^ ), one finds that Grover's output state vector after time 
T lays in a 2-dimensional sub-space, spanned by the Hadamard operation on the initial state 
of zeros W\0) and the marked state |m) (Note that rotates the state vector |0) in the 
2-dimensional sub-space spanned by |0) and iy|m) while = / is the identity operator). 

The unitary noise in the Hadamard transforms add perturbation elements to Grover's 
operator Go{T) [see Eq. (I28j) ]. These elements consist of the noise generators At and Ct, 
powers of the Grover iteration Q and selective inversion Jq. Executions of Q and Jq do 
not remain confined to a certain 3-dimensional sub-space. However, any operation of either 
At or Ct on a basis state vector of the form \y) or {\y) = \yn-i,yn-2, ■ ■ - yo) where 

yk = 0,1), extends the superposition by n additional independent state vectors. The higher 
the order / of the perturbation component Gi{T) [see Eq. (I25|) ]. the larger the number of 
executions of the noise generators. 

Consider a realization of the noise in which all the a's and the 7's are of order e. Table H] 
shows the vectors which appear in the superpositions produced by the leading perturbation 

15 



components of Grover's expansion ()28p with their corresponding order of the noise. An 
optimal time T = Tq, in which the noiseless Grover's operator returns exactly the marked 
state |m) is assumed (i.e. G'o(To)|0) = |m)). The vectors and {\Oik)}o<k<n-i 

are the n vectors that are respectively different from the marked state \m) and the state of 
zeros |0) in their k'th bit only, the vectors {\rh2k,k2)}o<k^<k2<n-i and {\02k^k2)}o<k^<k2<n~i 
denote the C2 vectors that are respectively different from those vectors in the /ci'th and 
fc2'th bits and so on ( C" is the binomial coefficient CJ^ = jff^^^zfy )- 

Given the marked state |m) we define a neighborhood class / which includes all com- 
putational basis states \x) that are different from the marked state |m) in exactly / bits, 
e.g. consider a 5 qubit system with marked state |m) = |17) = 1 10001), one finds the state 
\x) = |31) = |11111) in the 3'rd neighborhood class. 

Table lU clearly shows that in spite of T being the optimal time Tq, the noise causes a flow 
of probability from the marked state |m). (In case that T is not optimal, the probability 
flows out of the original 2-dimensional Grover's space, spanned by \rh) and iy|0).) 

A certain portion of the probability "diffuses" from the marked state \rh) to its neighbors 
{|"^ifc)}o<fe<„_n {|"^2fcifc2)}o<fci<fc2<n-i ^^^-^ where any perturbation component of order / 
(G'i(T)) contributes amplitudes of order O(e') to all the states within neighborhood classes 
of order I' < I. 

Another portion of probability involves the Hadamard operator and hence flows uniformly 
to all computational basis states. Since the first component which includes a Hadamard 
operator is O(e^), the contribution of this flow to each computational basis state is 0{^). 

Clearly in the limit of large A^, the "diffusive" flow into states of small / neighborhood 
class (of order 0(e')) is much stronger than the uniform flow (of order 0{^^)). We call these 
states |m)'s near states. All other states are considered far states. There is a trade-off; on 
one hand the probability to measure a certain near state is larger than the probability to 
measure a certain far one. On the other hand there are much more far states {0{N) far 
states where A^ = 2") than near states {0{n'') near states, where I is small). 

We will now quantify the performance of Grover's search algorithm in the presence of 
unitary noise. To this end we will calculate the average probabilities to measure certain 
basis states \x) = x„_2, . . . Xq) at the optimal measurement time. Without noise, 

at the optimal measurement time the probability to measure the marked state is po = 1- 
Thus, the probability of any other state is zero. According to Eq. (jHl) using the matrix 
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representation in Eq. (jHO)), the unperturbed probability of measuring an unmarked basis 
state (i.e. \x) ^ |rn)) at a certain time T is: 

Po(x,T) = |(a:|Go(T)|0)|' = |^ cos(^T) + 0(i)|' (56) 



where = 2" is the total number of computational basis states. Therefore the optimal 
measurement time is given by: 

Adding a perturbation to Grover's quantum search operator, a second order expansion 
of the probability to measure a state \x) at time T gives: 

p(x,T) = |(x|G(T)|0)|' = |(x|Go(T)|0)|' + AKx,T) (58) 



where: 



Ap(a;,T) = R{x,T) + |(a;|Gi(T)|0)|' (59) 



and 

R{x,T) = 2{x\Go{T)\0)Re{{x\Gi{T)\0) + {x\G2{T)\0)}. (60) 

The normalization condition dictates that the probability to measure the marked state \m) 
at any time T satisfies p{rh,T) = 1 — '^x=/=fnP{x,T). When T is around the optimal time 
To, the probability to measure the marked state approaches unity. Therefore, all the terms 
that are much less than 0{jj) and appear in the — 1 probabilities p{x,T) to measure 
unmarked states \x) respectively, become negligible. 

For time T close to Tq an expression of the form T = Tq + t where \t\ = 0(1) ^ To = 
0{^/N) can be written. Hence, for any unmarked state the estimation: 

|(x|Go(T)|0)| = -^1 cos(^T)| = -^1 sin(^t)| = ^\t\ = O (^1) (61) 

— 2 

can be made. Therefore, |(a;|Go(T)|0)| = 0{^) ^ 0{j^) is negligible. 

On the other hand, a small perturbation around the original Grover operator Go{T) is 
assumed. Namely, |(x|G'2(T)|0)| < |(x|Gi(T)|0)| < 0(1), so that |i?(x,T)| < O (^) and 
is therefore negligible. 
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Hence, the mean probability of measuring a computational basis state \x) at time T 
around the optimal time Tq is given by: 

{p{x,T)) ^ ^ . _ 2 (62) 

i l-E.^^((l(^|G'i(T)|0)| )) ifx = m 

where |m) is the marked state, () denotes the averaging on the noise and 

T 

(^|Gi(r)|0) =z{J](a;|Tyg^-*(A+i + Joa/o)Q*|0)} (63) 

t=o 

using the definition of Ct=o = 0. 

In the Appendix, we apply the matrix representations which appear in Eqs. (j5(Jj) and 
(jH^ as well as the noise matrix elements in Eqs. (fT!Hl - (f^ and (j22I)-(I2ni) to calculate the 
leading terms of the mean measurement probabilities at the optimal measurement time in 
the limit of weak noise and large A^. We assume that the two Hadamard transforms, which 
appear in each Grover iteration, are implemented with similar hardware. We also assume 
that the noise is unbiased, or alternately we refer to the case in which the exact values 
of the bias elements are known (a-priory) and can be shifted away in every iteration of 
the Hadamard gate. For further simplification we only consider an isotropic noise. Thus, 
statistical moments of the noise defined in Eqs. (jlSp and (fTB|) takes the form: 

{a,) = (7m) = (64) 

and 

(Sa^Sa^) = {S'j^S'j^) = 5^^e^ (65) 

where 6^^, is the Kronecker's delta function (yU, = 1,2,3) and e is the isotropic noise's 
standard deviation. It is shown in the Appendix (Eqs. ()A38|) - ()A40|) ) that for any given 
distribution of the noise, the averaged probabilities at optimal measurement time are ap- 
proximated in the limit of large A^ and small e by: 

Po = 1 - -TmVNe^ + 0{n^Ne^) (66) 



Pfar = ^nnVNe^ + 0{n^Ne^) 



(67) 
(68) 
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Po is the mean measurement probability of the marked state \fh). Pi is the averaged prob- 
abihty to measure a near computational basis state which lays in the first marked state's 
neighborhood class (and is hence different from the marked state in a single bit only). Pfar 
is the averaged probability to measure a far computational basis state that differs from the 
marked state in more than one bit. 

This result reveals some flexibility which enables us to find the marked state even when a 
searching error has occurred. In case of weak noise (i.e. if e is small), the mean probabilities 
to measure near and far states are of the same order. With an approximate probability of 
P^^P^ ~ I' ^ searching error yields a near state which is different from the marked state 
in a single bit only. By flipping the bits of the measured near state one at a time, the 
marked state can be reconstructed after at most n steps. Similar behavior also exists in the 
non-isotropic case, although the mean probabilities no longer have a compact form. 

A re-scaled standard deviation magnitude: 



may be considered as the parameter of the problem. The limit of 77 ^ 1 (e <^ J- — ) implies 
for weak noise in which the above approximations are valid. The limit of 77 ^ 1 (e ^ — ) 
is the strong noise limit where the noise completely destroys the quantum search. In this 
case the measured state is randomly taken from the total number of N possible measured 
states, so that Pq ^ and Pi — > In between those extremes we refer to the noise as 
moderate. 

In the next section we present numerical results which support the analytical predictions. 
We evaluate the performance of the unbiased and isotropic noisy quantum search also in 
case of moderate noise, where the limit of weak noise is not valid anymore. 

IV. THE PERFORMANCE OF THE NOISY SEARCH ALGORITHM - NUMER- 
ICAL RESULTS 

We first show simulation results confirming the predictions of the analytic approximation. 
The simulations also enable us to study the effect of larger values of e (the standard deviation 
of the noise) . These results give rise to new search strategies which we discuss in section V. 

The simulations were written in C-|— I- and FORTRAN, on 1686 machines running Red Hat 




(69) 



19 



Linux. We applied isotropic unbiased Gaussian unitary noise by transforming the output 
of a uniform random number generator to be Gaussly distributed. The random Gaussian 
variables ai, and 7^ (/i = 1,2,3) are the coefficients of the Pauli matrices as defined in 
Eq. (jl4j) . Note that although these coefficient are limited to the range between —it and 
vr, the standard deviation of the noise satisfies e ^ 1, thus the corrections to the Gaussian 
distribution are negligible. In order for the results to be statistically sound they are averaged 
over sufficient number of runs. The statistical error is estimated by comparing the results 
of different runs with identical parameters. 

Fig. HI shows the dependence of Po,Pi and Pjar on the standard deviation of the noise e, 
for relatively small values of e. The predictions of Eqs. (j67|) and (|68|) . are compared to 

the simulated results. The prediction is valid for small values of e (e < 0.006 ^ 0.166.^^^—^, 
for number of qubits in the register n = 12). In the presence of noise with larger standard 
deviation, the next terms in the approximation are no longer negligible, and the simulated 
results start to deviate from the predicted ones. 

Fig. El shows the same dependence for larger values of e. Results are plotted for system 
sizes of 8 and 15 bits. Note that a re-scaled standard deviation 77 = a/ n^/Ne axis is used 
as predicted in Eq. ()69|) . Indeed, this scaling makes the two graphs diverge. This further 
confirms the validity of the approximation. 

We can roughly divide the noise deviation level into three regions: 



Weak noise deviations: r] < 0.166 (e < 0.166^ / —j^)- This level of noise was previously 



rWN' 

discussed. 



Strong noise deviations: ri > 1.4 (e > 1.4w^=). At this level of the noise Grover's 
algorithm should completely fail and the probability distribution among the basis 
states should be uniform. This means that for large e, Pq ^ , and Pi ^ ^ which is 
verified by the simulation. 



Moderate noise deviation: 0.166 < < 1.4 (0.166^^ < e < 1.4^^). The exact 
behavior of Grover's algorithm under this level of noise and the specific boundaries of 
the moderate noise region are not predicted by the analytic approximation. However, 
the simulation enables us to study this C3iSG clS well. As the noise increases, probability 
"diffuses" from the marked state to the near states and Pi increases. Gradually, 
the standard deviation of the noise becomes too large, which makes the probability 
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"diffuse" to all states, near and far, uniformly. In this case the vast number {0{N)) 
of far states overcomes the few {n) near states, and Pfar approaches one. 

Note that there is no noise level for which Pq is negligible while Pi is not. This is a results 
of the effect of high order terms in the approximation. Recall equations (jUUI) and ()67|) . For 
Pi to be of the order of 1 we need = O (^^^T^j) and in this case the higher order terms 
are 0(1) as well, which implies that Pq is not necessarily negligible. 

To demonstrate the flow of probability from the marked state to near and far states, we 
divide the computational basis vectors into classes of neighborhood to the marked state. 
Given the marked state vector \m) we have defined its Uth order neighborhood class as the 
set of C" = nf^^^iy^ state vectors that differ from |m) in exactly / bits. Let Pi, I = 0, . . . ,n 
denote the probability to measure any of the state vectors of the Uth neighborhood class. 
Fig. Elshows the distribution of probability among the P;'s for different values of rj. For small 
1], the marked state (neighborhood class of order 0) should be measured with probability 
1. The remaining neighborhood classes should have zero measurement probability. For 
large r], the distribution of probability should be uniform among basis state vectors, which 
should reflect in a binomial distribution of probability among the neighborhood classes 
(proportional to their size). The transformation between the two extremes is determined by 
the two probability flow types caused by the noise: the "diffusive" and the uniform. First Pi 
increases due to the "diffusive" part of the flow. Then, the middle (around |) neighborhood 
classes become noticeable due to the number of elements in these classes. Yet, there are 
cases (0.803 < rj < 1.069) where Pq is quite small, but Pi and P2 are still relatively large. For 
even larger values of 77 we have a semi-binomial distribution distorted toward the low order 
neighborhood classes, which eventually (for strong noise) becomes an almost pure binomial 
distribution. 



V. ALTERNATIVE SEARCH STRATEGIES 



One can identify two search strategies using classical or quantum search: 

• Classical search until marked element is found. This strategy has an average run time 
of 

'^classical ('^^) 
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steps. 

• Quantum search followed by a single classical verification step until marked element is 
found. This is a geometric procedure with success probability Pq and has an average 
run time of ^ 

-TO 

where a single classical computation step is performed in one time unit and a single 
quantum (Grover) iteration is performed in Tq time units. 

The property of "diffusion" of probability from the marked state to the far states through 
the near ones raises the possibility of using hybrid (quantum and classical) search strategies. 
These strategies are based on classically searching the marked element starting from the 
state measured after applying Grover's algorithm, and going over the states according to 
the class of neighborhood to which they belong. The general hybrid strategy has a parameter 
I, and is defined as follows: 

1. Run Grover's algorithm, and measure the state of the register in the computational 
basis at the optimal measurement time {^\fNT^. 

2. Repeat for j = 0. . .i: Classically search the marked state among the states of the 
j'th-order neighborhood class until it is found. 

3. If the marked element was not found, go back to step 1. 

Naturally, the average number of quantum and classical operations required in order to find 
the marked element depends on the noise (which controls the probabilities P^), and on the 
value of the parameter I. The effectiveness of each strategy is also governed by the time it 
takes to perform a single quantum computation step, Tq (which is the time required for a 
single Grover iteration). 

In order to analyze the performance of each strategy we define (7]) to be the average 
time required to complete an entire search of the marked element using the hybrid strategy 
with parameter /. We also define Tli to be the probability that the marked element is found 
in a search of the first / neighborhood classes. Obviously, 
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Let Ti denote the average time required to complete a single execution of steps 1 and 2 
above in case that the marked element is found, and Ti denote the time required to perform 
a single execution of these steps in case the mark element is not found. These are given by: 



I 

TT 



Ti = -VNT, + J2cr (73) 
i=o 



where C" is the binomial coefficient C" = jfj^iiy , and 



Poi^VNr, + 1) + Pii^VNr, + 1 + + PslyV^r, + 1 + + ^-C^) + ... 



' j=i \i=i / 



(74) 



This is the calculation of the mean time, normalized by 11;. 
We can now express (T/) in terms of Ti,Ti and 11; : 

(Ti) = UiTi + Ui{l-Ui){Ti + fi) + Ui{l-Ui)\Ti + 2fi) + ... 

oo oo 

= UiTi ^(1 - UiY + urn J2 Ki - ^lY 

r=0 r=Q 

= Ti + ^^Tj. (75) 

The average time of the optimal strategy is given by 

(T)=m^in((T,)). (76) 

Of course we must have (T) < y. Otherwise the quantum search is completely inefficient, 
and the regular classical search yields better results. The best strategy parameter lopt is: 

I opt = arg ^min (T/) j , (77) 

where arg denotes the argument of the minimum function, i.e. the value of / for which (T;) 
is minimal. 

We work under the (strict) assumption that Tg = 1, which is that a quantum computation 
step requires the same amount of time as a classical computation step. Fig. 0] shows 
the optimal strategy averaged time (T), calculated from the data presented in Fig. 01 
as a function of the noise's re-scaled standard deviation t] = a/ ny/Ne It also shows the 
performance of the two trivial strategies mentioned above (Eqs. 1701 and I7T| . It is apparent 
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that the required time is decreased by up to a factor of 7. It is especially important to 
note that for noise deviation rj > 1.4 (e > 1.4 J- — =) the original Grover quantum search 
is useless. (At the point t] ^ 1.4 the curve of the naive quantum searching time intersects 
with the line which is the averaged time of a classical search). On the other hand, 
the optimal hybrid strategy can be used to substantially decrease the required time up to 
T] ^ 1.9 (e ^ 1.9 J- — =). It is obvious that under large noise deviations the optimal hybrid 
strategy requires ^^^^Tq + ^ — y to complete the search (this is the time required for 
an initial quantum search followed by a classical search of the entire search domain). As it 
approaches this limit ( ?7 > 1.9), the dependency of the time required by the hybrid strategy 
becomes a concave function of rj. Hence, it is for this noise levels only ( 1.4 < < 1.9 ) that 
the hybrid strategy yields significant improvements. Above this level of noise, the classical 
search is at least as efficient (and does not require the initial quantum search). 

Table mi shows the chosen strategy lopt calculated for each level of noise. It can be seen 
that as the noise's standard deviation increases, the chosen / increases as well. This means 
that it is worthwhile to classically search increasingly more distant neighborhood classes for 
the marked element. For large t], lopt = 20, which means that the chosen strategy is not 
better than the naive classical search. 

This analysis shows that the property of "diffusion" of probability from the marked state 
to the near states can be used to obtain a hybrid search strategy which is more efficient than 
the naive quantum search strategy. Note that the results above were calculated for Tq = 1. 
If we make a more reasonable assumption Tq ^ 1 (i.e., that a single quantum computation 
step takes more time than a classical one), the improvement would be considerably better. 
In addition, the hybrid strategy enables an efficient quantum search under noise levels for 
which the naive quantum strategy fails. 

VI. SUMMARY 

We have studied the effect of unitary noise on the performance of Grover's quantum search 
algorithm. This type of noise may result from tiny fluctuations and drift in the parameters of 
the (quantum) components performing the computation. The resulting operations are still 
unitary, but not precisely those assumed in the design of the algorithm. In the analysis we 
focused on the effect of an unbiased unitary noise in the Hadamard gate W. For simplicity 
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we further assumed that the noise is isotropic as well. The Hadamard gate is an essential 
component in each iteration of the quantum search process. The gate W was replaced by 
a noisy Hadamard gate U, whose parameters are distributed around those of W according 
to an unbiased, symmetric probability distribution. The noise level was characterized by its 
standard deviation. 

It was found that for noise levels greater than 0{n~^ N^^), Grover's algorithm becomes 
inefficient. (Here n denotes the number of qubits in the register and = 2" is the number 
of computational basis states). The nature of the flow of probability out of the marked 
state, which is caused by the noise was also investigated, analytically and numerically. A 
phenomenon of "diffusive" flow of probability to the marked state's near neighbors, that are 
different from the marked state in only a few bits, was observed. This feature of "diffusion" 
gives rise to new hybrid search strategies which are shown to improve the effectiveness of 
the entire search procedure, both in terms of the time required to complete the search, and 
in the noise level under which a quantum search is still more effective than the classical 
search. The use of hybrid strategies was found efficient even under noise levels for which 
naive re-execution of Grover's quantum search (until the marked state is found) fails. 
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APPENDIX A: CALCULATION OF MEAN PROBABILITIES 

We are aimed to obtain an expression for the averaged probabilities to measure a state 
\x) around the optimal measurement time in the presence of one qubit unitary noise in the 
Hadamard gate. We begin with Eq. (jU^ . As one can see (x|G'i(T) |0) is built as sum over 
time index t of matrix elements produced by 4 consequent operations on the initial state of 
zeros |0): Q*, (^t+i + loCtIo), <5"^~* and W in that order. Using the matrix representations 
of and Jo in Eqs. fl5()|l and ()54|1 . as well as the Hermitian noise generators presented in 
Eq. (jSSI), the following operations can be expanded: 




(Al) 
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n-1 



+ /oC'i/o)|0) = (0|(A+i + Ct)\0)\0) + J](Oifc|(A+i - ^010)101^) (A2) 



k=0 



(it+i + ioCJo)W\m) = (m|iy(A+i + Ct)W\m)W\m) 



Y,{mik\W{A+, + Ct)W\m)W\m,k) + Oq (A3) 



k=0 



Q^-'\0) = cos(cj[r - t])\0) + sm(cu[T - t])W\Th) + 6 (A4) 
Q^-'\Oik) = i-lf'%k) + O (A5) 
Q^-'W\m) = ~sm{uj[T -t])\0) +cos{uj[T -t])W\m) + (A6) 

Q^-'W\m,k) = {-lf-'W\m,k) + 6 (A7) 

where Oifc and fhij, are the binary strings which are respectively different from the ini- 
tial string of zeros and the marked string fh only in their fc'th bit, = {^'^^ |0) + 
O (^"^j W\ffi), and Oq = O(-^)|0) where 0(e) denotes the order of the noise. Then, 
performing consequent substitutions of Eqs. ()A1|) - ()A7|) in Eq. ()63p while memorizing that 
(a;|iy|0) = and (x|W^|Oia:) = (—1)'^'''-^ (xfc is the binary value of the fc'th bit in the 
string X = x„_i, x„_2, . . . ,Xo) one finds that for any unmarked state \x): 

T 



{x\G,{T)\0) ^cos(a.t)cos(a;[T-t])(0|(A+i + a)|0) 

t=i 

T 

- J2 sin i^t) sin {uj[T - t]) {m\W {At+i + Ct)W\m) 
t=i 

T n-1 

+ 5^cos(.;t)(-l)^-*[5^(-l)-'=(0i.|(i,+i-(7,)|0)]} 

t=l A:=0 
n-1 T 

+z {^(x|mifc)^sin(^t)(-l)^-*(mifc|iy(it+i + a)W^|m)} 
+0 



fc=0 t=i 
'Tne 



N 



-^j+g(x|mi,)0(rne), (A8) 
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where according to Eqs. (fT^ -(P ^ : 

n-l 

(0|(^+i + Ct)\0) = 5^{aofc(t + 1) + a^kit + 1) + 7ofc(t) + 73fc(t)} (A9) 



k=0 



n-l 

(m\W{A,,+i + Ct)W\m) = Y^{aok{t+l) + {-ir'a^kit+l) +lok{t) + (-l)™'=73,(t)} (AlO) 

fc=0 



(Oifc|(it+i - Ct)\0) = aik{t + 1) + ia2k{t + 1) - 7ifc(t + 1) - ^l2k{t + 1) (All) 



(mifc|W^(A+i + Ct)W\m) = ask{t + 1) - ^(-l)"^'=«2fc(t + 1) + 73fc(t) - 2(-l)'"'=72fc(t). (A12) 

We immediately observe that the order of {x\Gi{T)\0) depends on the Hamming distance of 
the state \x) to the marked state |m). The order of "far" states, which denote binary strings 
different from the marked string in more than one bit is 0{\/N) times smaller than the 
order of "near" states whose binary strings are different from the marked string in exactly 
one bit. 

" — 2 

According to Eq. we focus on computing the mean value of |(2;|Gi(T)|0)| around 

the optimal measurement time. For this purpose we evaluate trigonometric sums which 
appear in (|(a;|G'i(T)|0)| ), assuming that uT = | + O where u ^ 

For any small angular frequency uj and arbitrary phase ip: 

El J. \i -,\T-t cos(co'r + (^) + cos(to'[T + 1] (-1) {cos{u + ^p) + COS cp} 
cosM + ^)(-l) = 2(1 + COS..) 

= 0(1) (A13) 



t=l 



Therefore: 



T 



J2 cos{ut){-lf-' = 0(1) (A14) 



t=i 

T 



sin(cut)(-l)'^"* = 0(1) (A15) 



t=i 



and: 



^cos2(cut)cos(cu[r-t])(-l)^-* = 0(1) (A16) 



t=i 
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sm{u!t) cos{ujt) cos{u[T - t]){-lf-' = 0(1) (A17) 



t=i 

T 



sin(a;[r-i])(-lf~* = 0(l). (A18) 



t=i 



(Note that any multiplication of sines and cosines can be reduced to sums of sines and 
cosines) . 

On the other hand: 

Ysin^ut) - - - + 1{1 _ cos{2uT)} = |V]V + 0(1) (A19) 



t=l 



y cos2(u;t) = - + _ 1(1 _ cos(2a;r)} = ^V^V + 0(1) (A20) 



Vcos(u;i)cos(u;[r-i]) = - cos(u;r) + ^^^^^^^ = lv]v + 0(l) (A21) 
^ ^ > ^ ^ ^' 'I ^ > 2 tana; 4 ^ ^ ^ ' 

T sin(a;r) 1 



T 

^sin(cut) sin(cu[r - t]) = ~ cos(cuT) + ^ii!^ ^ ^v^+ (A22) 

t=i 



ycos''{u;t)cos\u;[T-t]) = -{1 + 1 cos (2a;r)} + ^^^^^^{ + ^} 

^ ^ ^ ^ ^ 4I 2 ^ 4 ^tan(2a;) tana;'' 

= ^V^+0(1) (A23) 



ysin2(cut)sin2(^[r-t]) = ^{l + icos(2^r)} + ^^^^^^{^— -— ^} 

4^2^^-' 4 Han(2u;) tano;^ 

= ^^ + 0(1) (A24) 

E'^ / rm / rm ^ / m\ sin(2a;T) 
cos(cc;t) cos(a;[T — tj) sm(ci;t) sm(a;[T — tj) = — — - cos(2ci;T) + — — 
8 8tan(2cx;j 

= ^v^+0(l). (A25) 

We have now completed all the needed expansions for expressing the mean probability to 
measure a certain unmarked state \x) around the optimal measurement time. By multiplying 
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Eq. ()A8|) with its complex conjugate and averaging the value of |(x|(j'i(T)|0)| using the 
first and the second moments of the noise (see Eqs. (fT3j) - (fTBj) ) while leaving the leading 
terms only, one obtains that the mean probability to measure an unmarked state \x) around 
the optimal measurement time is: 

p, ^ l^{{Sa2') + {6a^) + (572') + (573')} + 

o 

^^[(«3) + (73)]-/(m)[(ai) + (7i)]}V 

((«3) + (73))' + ((«2) + (72))' (A26) 

in case that is a first order neighbor of the marked state \rh) (i.e = Imi^), is different 
from the marked state only in its /c'th bit) and 

Pfar {^[(a3) + (73)]-/(m)[(ai) + (7i)]}' + 

f{m)[{6a,6as) + {6jiSrs)]} (A27) 

in case that \x) is far state (i.e. the state \x) differs from the marked state |m) in more 
than one bit). The a^'s and 7^'s {fi = 1, 2, 3) denote the real stochastic variables taken from 
any arbitrary distribution which characterize the noise, (The moments of and 70 vanish 
because they act as a global phase), n is the number of qubits in the register, /(m) is the 
difference between the number of zeros and the number of ones that appear in the binary 
representation of the marked string m and N = 2^. 

In the noiseless search algorithm, the n-qubit Hadamard transform (that consists of n one 
qubit gates) appears twice in each iteration. These n operations are expected to be imple- 
mented with similar hardware. The noise characteristics of all one qubit gates are expected 
to be similar but with no correlations between each other. The noisy unitary operators Ut 
and Vt at any Grover's iteration t can therefore be considered as direct multiplications of n 
one-qubit unitary and stochastic operators of the same statistical behavior. Thus, according 
to Eq. (fT^ the statistical properties of Ck{t) and Wkaic{t)wk are alike for any time index t 
and qubit index k. Here Wk denotes a single-qubit operator, acting on the fc'th qubit i.e.: 

Wk = + ^3fc) 

where aik and ask are Pauli operators which act on the /c'th qubit. Using the expansions of 
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ak{t) and Ck{t) by Pauli operators (see Eq. ()13|)) with the aid of the identity; 

{vl ■ (y){v2 ■ (y) = {vi ■ V2)cro + ^(^1 x ?}^) ■ a 

where vl and V2 are 3-dimensional vectors of real numbers, a = ((Ti, (3"2, (T3) is the vector of 
Pauh operators, o"o is the identity operator, and ■ and x denote scalar and vector products 
respectively, one obtains the following relations between statistical moments: 

(71) = ("3) (72) = -(«2) (73) = («i) (A28) 

and 



(571572) = -{5a25a3) (571^73) = (5«i5«3) (572^73) = -{5ai5a2). (A29) 

Moreover, further simplification can be done, if we assume that the noise is unbiased, or 
alternately that the exact value of the bias of the noise is known and can be shifted in every 
operation of the Hadamard gate. We also assume that the physical system which realizes 
the quantum computer is happens to be isotropic. Thus, the statistical moments of the noise 
become: 

("m) = (7m) = 

and 

where S^^, is the Kronecker's delta function (yU, = 1,2,3) and e is the isotropic noise's 
standard deviation. Then, substituting these relations in Eqs. ()A26|) and ()A27|) one gets: 

Pi ^ -VNe^ (A30) 
2 



and 



Since there are n first order neighbors to the marked state |m) and N — n — 1 far states, we 
find that: 

Pi = npi ^ ^n\fNe^ (A32) 
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is the mean probability to measure a first order marked state's neighbor in the optimal 
measurement time and 

Pfar ~ Npf^r ~ ^vrnViVe^ (A33) 

o 

is mean measurement probability of an unmarked far state (in the optimal time as well). 
The mean probability to measure the marked state is given by the normalization condition: 

Po = l-Pi-Pfar^l- lunVNe''. (A34) 

o 

In order to evaluate the limits in which the above approximations are still valid, we have 
to estimate the order of the residual component produced by higher terms of the noise. In 
our calculation we have focused on the case of an unbiased noise. Therefore, only even 
powers of the noise standard deviation e appear in the mean probabilities expansions, so the 
next term in the above approximations is 0{e'^). By taking those terms into account, Eq. 
fl62|) has the form: 

(p(x,T)) ^ (|(a;|Gi(T)|0)|V (1(^162^10)1') + 

(((x|Gi(T)|0)(x|G3(T)|0)*)) + (((x|Gi(T)|0)*(x|G3(T)|0))) (A35) 

where {p{x,T)) is the mean probability to measure a certain computational basis state at 
time T around the optimal measurement time, and Gi(T), 62 (T) and 63 (T) are Grover's 
perturbation components of the first, second and third order respectively, as defined in Eq. 

An explicit calculation shows that in case that the measured state \x) is a marked state 
first order neighbor (i.e. \x) = \rhik) is different from the marked state \m) in the fc'th bit 
only): 

O{{x\G^{T)\0){x\Gs{T)m = O(|(x|G2(T)|0)|') = 0{TV) (A36) 
while in case that is far from the marked state \rh): 

2 T^n^e^ 

O{{x\G,{T)\0){x\G,iT)m = O(|(x|G2(T)|0)| ) = 0(^^). (A37) 

Thus, due to the fact that there are n first order neighbors of the marked state and N—n—1 ^ 
far states, one finds that the mean measurement probabilities at time T around the 
optimal measurement time (such that 0{T) = 0{\/N)) are given by: 

Po = 1 - InnVNe^ + 0{n^Ne^) (A38) 
8 
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^1 = 7T 




(A39) 



5 



(A40) 



ar 



8 



where Pq, Pi and Pfar are the mean probabihties to measure the marked state, a first order 
neighbor of the marked state and a far state respectively. 
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TABLE I: The list of vectors which appear in the superpositions produced by the leading pertur- 
bation components of Grover's expansion H28|) with their corresponding order of the noise. The 
first row of the table denotes the noiseless Grover's search result at the optimal time Tq, which 
satisfies by definition Go(To)|0) = \m). The second and third the row lists the vectors which appear 
in the superpositions produced by the first and the second order perturbation components Gi and 
G2 respectively. 



Perturbation component 


Order of the noise 


Vectors in the superposition 




0(6°) 




Glim) 




{|?Tllfc)}o<fe<„_i 
H^|0>, {W\Olk)}o<k<n-l 


G2{Tm 




N>> {I^U-)}o<fe<n-l' {l'^2fcifc2)}o<fei<fc2<n-l 
W\0), {W\0^k)}o<k<n-V {W\02k,k,)}o<k,<k,<n-l 









°In case that T is not the optimal measurement time the vector 14^ |0) also appears in the list. 

TABLE IL The optimal strategy lopt for different levels of noise deviation. The calculations are 
based on the data shown on Fig. 13 (number of qubits is 20). The noise re-scaled standard deviation 
is ?7 = yru/Ne. 



V 


0.0053 


0.138 


0.271 


0.404 


0.537 


0.670 


0.803 


0.936 


1.069 


^opt 





1 


1 


1 


1 


2 


2 


2 


2 


V 


1.202 


1.335 


1.468 


1.607 


1.734 


1.867 


2.001 


2.133 




^opt 


2 


2 


2 


3 


3 


3 


19 


20 


20 
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